Dynamical control and novel quantum phases in impurity doped linear ion crystals 
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We explore the behavior of the phonon number distribution in an heterogeneous linear ion crystal. 
The presence of ion species with different masses changes dramatically the transverse energy spec- 
trum, in such a way that two eigenfrequencies become non-analytic functions of the mass ratio in the 
form of a sharp cusp. This non-analyticity induces a quantum phase transition between condensed 
and conducting phase of the transverse local phonons. In order to continuously vary the mass ratio 
we adiabatically modify a locally applied laser field, exerting optical dipole forces which reduces the 
effective mass. 

PACS numbers: 03.67.Ac, 37.10. Ty, 64.70. Tg, 73.43.Nq 



I. INTRODUCTION 

Trapped ions are one of the most attractive physical 
systems for implementing quantum computation |l| and 
quantum simulation Q. The ability to control and mea- 
sure the internal and external degrees of freedom with 
high accuracy allows for experimental implementations 
of various quantum gates [3H6[ and quantum comput- 
ing protocols • Trapped ions represent a convenient 
system for the simulation of many-body effects, such as 
quantum phase transitions in spin systems [lol [Tl|. atoms 



in optical lattices [12[ and coupled cavity arrays [13[ . Re- 
cently a simulation of the Dirac equation for a free spin- 
1/2 particle [3| has been successfully performed with a 
single trapped ion [TBj . 

While small, homogeneous ion crystals have been ex- 
tensively studied, the novelty of this paper is elucidating 
the richness of quantum phases in ion crystals doped with 
a second ion species. Depending on the mass ratio, the 
crystal allows for the observation of two quantum phases 
of the transverse local phonon number distribution. For 
a heavy impurity ion all phonons are condensed by the 
latter and the variance of the transverse local phonon dis- 
tribution is zero. At the same time zero average phonons 
with vanishing variance and correlation are observed at 
all other sites. In the following this state is referred to as 
the quantum condensed phase. In the opposite case, if 
a lighter impurity ion is added, the system is in a quan- 
tum conducting phase [12j where the phonons are redis- 
tributed among all ions showing non-zero variance and 
correlation, while the average phonons of the impurity 
ion is significantly reduced. There are two ways to vary 
the mass ratio in heterogeneous ion crystal. The first one 
is by doping a mono-species ion crystal with one impu- 
rity ion, a technique which has been applied for ion fre- 
quency standards [l6l |. sympathetic cooling of molecular 
ions [171 ] , and for deterministic ion implantation [18[ . The 
second scheme overcomes this discrete, stepwise variation 
and uses optical dipole forces by non-resonant tightly fo- 
cused laser beams [19( (Fig. Q}. It is still advantageous 
to use heavier doping ions because the dipole force can 
only reduce the effective mass but to observe the phase 
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FIG. 1: (color online). Impurity-doped linear ion crystal con- 
sisting of ions with mass m and an impurity ion with mass 
M. The mass ratio p = M/m determines the two quantum 
phases. Applying a tightly focused laser field to the impurity 
ion in the transverse x direction makes it possible to continu- 
ously vary the effective mass ratio. Increasing the magnitude 
of the light field leads to a decrease of the effective mass. 



transition the critical point of unity mass ratio has to be 
crossed. Current ion trap technology allows to implement 
all crucial elements of the proposal: preparation of the 
initial state, laser-ion interaction driving the dynamics 
required for the quantum phase transition, and readout 
of the phonon distribution. 



II. TRANSVERSE ENERGY SPECTRUM 

In the following we first consider a harmonically con- 
fined ion crystal with N — 1 ions of mass m and one 
impurity ion of mass M at position jm without any fo- 
cused laser field. If the radial confinement is stronger 
than the axial one then ions are arranged in a linear ion 
crystal along the axial z axis and occupy equilibrium po- 
sitions z®. Since the axial trap potential is independent 
of the mass, the equilibrium positions zf = liii of the 
ions are independent of the composition of the ion crys- 
tal. Here I is the natural length scale and u.i are dimen- 
sionless equilibrium positions [20l ]. The transverse, dy- 
namic Paul confinement along x and y directions is gen- 
erated by an applied radio frequency (RF) quadrupolc 
field. In the following we only consider one transverse 
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direction x. The oscillation frequency in x direction is 
uj x = |e|L^riFCx / '(v^'tiOrf) , where £/rf is the amplitude 
of the RF voltage with frequency Orf, e is the electron 
charge and is a geometry factor. The transverse os- 
cillation frequency is additionally reduced by the axial 
trap frequency w z according to lo x — u> x (l — a 2 /2) 1 / 2 , 
in lowest-order approximation, where a = ui z /uj® [2lj |. 
For a small displacement the motional degrees of free- 
dom in x, y and z direction are decoupled. The normal 
mode eigenfrequencies ojk — u! x y/Xk and eigenvectors b k - 
are obtained from the diagonalization of a matrix 
with YliBijb* = Xkbj, where uj x is the transverse fre- 
quency of a single ion with mass to. The real and sym- 
metric iV-dimensional matrix Bij is determined by the 
harmonic expansion of the external trapping potential 
and the Coulomb interaction between the ions; it reads 

M 



Bi 
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where = M/m denotes the mass ratio of both ion 
species. By expressing the position and momentum op- 
erators in terms of normal phonon modes 

%« = -^X>L**' p^ = ^E fe L^ ( 2 ) 



the Hamiltonian takes the form H — J2k ^fc^fc + l/^)- 
Here X}~ = *JK/2nuJk~(j£ + is the normal mode posi- 
tion operator and Pk = i^Hrnu>k/2(a k — dfc) is the nor- 
mal mode momentum operator, while at and a& are the 
phonon creation and annihilation operators of the kth 
collective phonon mode and N k 
collective phonon number operator 

In order to visualize the transverse mode we plot in 
Fig.[2]the eigenfrequencies for a string of six ions, with 
an impurity ion at position ]m = 2 when the mass ra- 
tio is varied. As \x increases beyond unity, the frequency 
spacing between the lowest-lying energy mode (LL) with 
frequency ojn and all other modes increases, while for 
fx < 1 it tends towards degeneracy. In the transverse 
direction the highest energy mode, the center-of-mass 
mode (COM) with frequency uii, is nearly degenerate 
for fj, > 1. When fi decreases below the value /1 = 1 
the frequency spacing between COM mode and all other 
modes increases. In real experiments the mass ratio can 
only be changed to discrete values as shown in Fig. [3] 



a£dfc is the respective 
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FIG. 2: (color online). Normalized transverse eigenfrequen- 
cies ujk/^Jx = V%c when /i is varied. The highest en- 
ergy center-of-mass (COM) and the lowest-lying (LL) modes 
are shown. Different ion species are depicted as follows: 
/11 = 9 Be+/ 24 Mg+, M2 = 24 Mg+/ 40 Ca+, M3 = 40 Ca+/ 43 Ca+ 
and m = 40 Ca+/ 27 Al+. b) Decreasing a the LL and COM 
modes develop a cusp at fi = 1. 



However, later in the paper we show how to overcome 
this limitation by optical forces. 

As a is decreased the gap at the avoided level-crossing 
(/i = I) between the modes vanishes. Moreover, LL and 
COM modes develop a cusp at /x = 1. For sufficiently 
small a <C 1 the off-diagonal elements of B^ fi} can 
be neglected and the the eigenfrequencies uik behave as 
follows 



uii — > 



/' 



for fj, < 1, 




This implies that the COM frequency u)\ and the LL fre- 
quency ujn, respectively, have discontinuous derivatives 
at u = 1. Therefore, the non-analyticity of the ground 
state energy indicates a quantum phase transition [23| . 
We will exemplify the properties of the two quantum 
phases by calculating the phonon number distribution, its 
variance and correlation. These quantities are all experi- 
mentally accessible by laser-spectroscopy measurements. 



III. QUANTUM PHASES 

The characteristic feature of the quantum phase tran- 
sition is expressed clearly by the behavior of the local 
phonon number operator 



= 271 ' 

Mu° x (M) ^ 

2h 



PI 



2tvmuj x [ {m) 
s2 



P 



JM 



2hMuj° x (M) 



(3) 



Initially the mixed ion crystal is prepared by laser-ion in- 
teractions in the state with n phonons in the LL collective 
mode and the other modes are cooled to the ground state, 
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FIG. 3: Average phonon number distribution as described by 
Eqs. Q and (J7J. The impurity ion is placed in the center of a 
string of N — 15 ions, i.e. at position ju = 8. The quantum 
condensed phase is observed in a heterogeneous ion crystal 
with fourteen 40 Ca + ions and one 43 Ca + impurity ion, where 
n — 2 collective phonons are localized at the impurity site (a) . 
For fx — 40/43 corresponding to fourteen 40 Ca + ions and one 
43 Ca + impurity ion the system is in a quantum conducting 
phase with non zero phonon distribution (b). 

achieving the state \*S?) = 1 00 . . . n). We find that the av- 
erage local phonons per site (per ion) (fij) = (^f\hj\^f) is 
given by 

(4) 

First we will discuses the properties of Eq. ((4]) for \i = 1 . 
The results implies that in general Y) A n i) — n - The in- 
equality does not violate the conservation of energy since 
the sum of the local phonons energies is decreased due 
to Coulomb interaction between the ions. However, if 
a is kept sufficiently small (a <C 1) the transversal vi- 
brational modes are approximately degenerate ui^ s» cu® 
and the local phonons per site are (nj) — n(b^) 2 so that 
^2,j{fij) — n. In this regime where the Coulomb energy 
is much weaker than the trapping energy the linear ion 
crystal is in the quantum conducting (superfluid) phase 
and the phonon number n is conserved [12] . Now we 
will discuss the properties of Eq. Q in the two limits 
(i > 1 and fi < 1, respectively. As can be seen from 
Fig. [2b for fi > 1 and a <C 1 all modes are degenerate 
except ojn — w®/fj,. Moreover, we find that b^ M = 1 and 

bjf { N — using Eq. (JXJ) . This implies that the heavy ion 
is oscillating while the others are at rest. The average 
number of transversal phonons is given by 

(nj^ ju )=0, (n jM )=n. (5) 

Thus, for a heavy impurity ion, increasing the number of 
collective phonons in the LL mode leads to an increase of 
the average number of local phonons at the impurity site 
whereas all other ions remain in the vibrational ground 
state. Hence, the heterogeneous ion crystal is in a quan- 



FIG. 4: The variance of the local phonon number operator as 
described by Eq. ([7]). The lighter impurity ion is placed in 
the center of a string of ./V = 15 ions, i.e. at position jm — 8. 
For (i = 40/43 corresponding to fourteen 43 Ca + ions and one 
40 Ca + impurity ion the system is in a quantum conducting 
phase. The phase is characterized with non-zero variance for 
all sites except for the impurity ion where the variance is zero. 

turn condensed phase where any transversal quantized 
energy is condensed by the heavy ion. 

In the opposite limit of a lighter impurity ion we derive 
from Eq. {1} that bj M = 1, b^ 1 = 0. Because now only 
the LL modes are excited, the oscillations of the impurity 
ion are suppressed and the average number of phonons is 

(n^ jM )=n(bff, (h jM )=0. (6) 

Hence, the lighter impurity ion spreads the phonons to 
all other ions, such that the heterogeneous ion crystal is 
in a quantum conducting phase. Eq. ((5]) shows that the 
lighter impurity ion is in a vibrational ground state and 
is independent to the collective phonons n. The phonon 
distribution is shown in Fig. El where both phases are ex- 
emplified for the system of a mixed 43 Ca + / 40 Ca + crystal. 

The different quantum phases become evident also 
from the variance of the local phonon number operator 
Srij = ({fij))— {fij) 2 ) 1 / 2 . Again we will start with the case 

of /i = 1, where we obtain that Srij = ^J n (bf) 2 — n(b^) 4 . 

However, for the heavy impurity case (fi > 1) we derive 
Srij = for all sites. The vanishing phonon variance ap- 
pears to due the the localization of the collective phonons 
at the impunity site. If the impurity ion is lighter (jj, < 1) 
we obtain 

8n jftj „ = yfiOff ~ n(bf)\ 5n 3M = 0. (7) 

Now the variance increases with the amount n of phonons 
in the LL mode for all ions except for the impurity ion. 

Finally, we demonstrate the quantum phases by con- 
sidering correlations in the number of phonons Cy = 
{fiifij) — {hi){hj). For fi = 1 the ion crystal is in 
the quantum conducting phase with correlation Cy = 
—n(b^b^) 2 . Since in the quantum condensed phase 
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(/i > 1) the collective LL phonons are localized at the 
impurity ion whereas the remaining ions are in the vi- 
brational ground state we find that the correlation Cy 
for any i and j vanishes. In the quantum conducting 
phase (/i < 1) the lighter impurity ion is in the vibra- 
tional ground state and therefore the correlation Cy M 
remains zero. However, for all other ions the correlation 
is Cy = -n{bfbff. 

To summarize, the impurity doped ion crystal allows 
to distinguish two quantum phases depends of the mass 
ratio fi. In the quantum condensed phase (/i > 1) the 
collective LL phonons are localized at the impurity ion 
whereas all other ions remain in the vibrational ground 
state. The localization leads to vanishing of the phonon 
variance and correlation at all sites. In the quantum con- 
ducting phase (fi < 1) the lighter impurity ion is in the 
vibrational ground state and the collective LL phonons 
are distributed among the remaining ions. This phase 
is characterized with non-zero variance and correlation 
except for the impurity ion. 

In the following, we show how to employ an optical 
dipole force in order to change the effective mass of the 
impurity ion. In this way \x can be continuously swept 
through the quantum phase transition. 



IV. OPTICAL DIPOLE FORCE 

We assume that additionally to the trapping potential 
an optical dipole force is applied to the impurity ion in 
the transverse x direction. This force creates a harmonic 
potential Vd = Moj 2 x 2 m /2 with frequency oj s . The total 
potential is the sum of the two trapping potentials and 
the Coulomb interaction between the ions. The Hamil- 
tonian for the transverse ion motion is given by [221 ] 



m(w°) : 



N 

£ 



BjiXjXi (8) 



Here for simplicity we normalize the position operators 
Xj as Xj = £j and stj M = ^JJiXj. Using Eq. ([2]) the 
Hamiltonian can be expressed in terms of normal modes 
as follows 



H(t) = H - 



1 N 

2 £ [ mRk 



q XkX q 



Skq{PkX q — XkPq) 



k,q=l 



(9) 

The new eigenfrequencies uik and eigenvectors b^j are 



obtained by diagonalization of a matrix By = By + 
P 2 &ijM&jjM i with /? = w s /w°. Due to the time-dependent 
optical dipole interaction the normal modes and mo- 
menta are connected by the couplings Rk q — XyLi bj bj 

and Sk q — J2f=i bjbj- Here we denote with dot the time 
derivative, occurring because the magnitude of the dipole 
force is varied in time. 

The experimental scheme starts by preparing the crys- 
tal in state |\E') = |00...n). During an adiabatic in- 
crease of the trap potential created by the optical dipole 
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FIG. 5: (color online), (a) The transverse eigenfrequencies 
w/t for COM and LL modes as a function of the optical-dipole 
frequency oj s for a heterogeneous ion crystal consisting of five 
40 Ca + ions and one 43 Ca + impurity ion at the second po- 
sition. A quantum phase transition is observed by adiabat- 
ically sweeping cj s to higher values, (b) The local average 
phonon distribution for each ion (ru) — m (i = 1,2..., 6) 
as a function of cj s . As the magnitude of uj s is changed the 
mass ratio is reduced and reaches a critical point fj, e g — 1 at 
which the system undergoes a quantum phase transition from 
the quantum condensed phase to the quantum conducting 
phase. The quantum phase transition occurs approximately 
at lj s = 0.37oj2. 



force the normal modes become coupled. In order to 
suppress non-adiabatic transitions to the other modes 
we require the adiabatic condition to be fulfilled at any 
instant of time, \Sk q \ <C \wk(t) — uj q (t)\. Then the cou- 
plings Rkq and Skq can be neglected such that the Hamil- 
tonian ([9]) describes TV uncoupled harmonic oscillators 
with time-dependent frequencies Wfe. We assume that 
the frequency of the tightly focused laser field is vary- 
ing as uj s = Wj y/t. The energy spectrum of the linear 
ion string is now represented by the Demkov-Osherov 
model [24[, which describes the interaction between N— 1 
stationary states and one state varying linearly in time. 
This allows us to approximate the adiabatic condition as 
(Wj/w") 2 < \u)k — u q \ 2 T. This condition implies that the 
temporal energy change should be smaller than the split- 
ting to the next level to avoid transitions. Then the ion 
crystal remains in the initially prepared state |^) during 
the adiabatic change of the optical dipole force. To il- 
lustrate the effect of the optical force on the transverse 
spectrum we plot in Fig. [5^, the eigenfrequencies for 
the COM and LL modes as a function of lu s for five 40 Ca + 
ions and one 43 Ca + impurity ion at jm — 2. For suffi- 
ciently small a <C 1 the off-diagonal elements of By can 
be neglected and we obtain that the optical dipole force 
reduces the mass ratio as // e g = ((uj s /lo®) 2 + 
As the magnitude of /z e ff is changed, the system reaches 
a critical point /i = 1 for which the system undergoes a 
quantum phase transition from the quantum condensed 
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phase to the quantum conducting phase (see Fig. [5]d). 
Here we assume that the heterogeneous crystal is pre- 
pared in a state with n — 2 LL collective phonons. In 
the quantum condensed phase these phonons are local- 
ized by the impurity ion with no phonons at the other 
sites. For sufficiently small a the adiabatic condition re- 
quires a long interaction time and high value of u> s . To 
avoid this problem we could increase a to 0.1. Then 
the quantum phase transition occurs approximately at 
lu s = 2tt x 0.4 MHz and the adiabatic condition is ful- 
filled within 60 /is. The optical dipole trap can be created 
by laser detuning up to A = — 2ir x 300 GHz red of the 
resonance of the Si/ 2 *H> P1/2 transition of 43 Ca + . The 
frequency lj s can be achieved by laser power around 350 
mW with waist radius of wq — 5 /im. 

V. CONCLUSIONS 

In conclusion we have presented a realistic scheme for 
observation of phonon condensation and quantum phase 



transition of the transverse local phonons in a hetero- 
geneous ion crystal. Critical behavior appears when the 
effective mass ratio is changed continuously, which can be 
achieved by application of light-induced dipole force. The 
features of the phase transition are clearly visible even for 
a small number of ions, realizable with the current ion 
trap technology. Impurity-doped ion crystals offer the 
advantage that the laser light can be made to interact 
only with the impurity ion. Moreover, the local phonon 
mode readout of the impurity ion can be improved by us- 
ing laser frequencies selectively addressing its sidebands 
and internal electronic states [13j . The future extension 
of the proposed technique may use multiple impurity ions 
for heat, excitation and entanglement transport measure- 
ments, investigations of artificial Josephson junctions and 
separation of phases. 
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